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D ■ Abstract 

New developments in carlomat, a program for automatic computation of 
the lowest order cross sections, are presented. They include improvements of 
. the phase integration routines and implementation of extensions of the stan- 

dard model, such as scalar electrodynamics or the anomalous Wtb coupling 
\ including operators of dimension up to five. 

1 Introduction 

carlomat is a program for automatic computation of the lowest order cross 
^ sections, dedicated in particular for the description of multiparticle processes 

of the form 

P1+P2 -»> P3 + ■■■ +Pn (1) 

with the maximum number of external particles n — 12. In (OQ), particles have 
been symbolized by their four momenta in the centre of mass system (c.m.s.). 
The program is written in Fortran 90/95. It generates the matrix element for 
a user specified process together with different phase space parametrizations 
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which are used for the multichannel Monte Carlo integration of the lowest 
order cross sections and event generation. 

Version 1.0 of carlomat was released 2 years ago pp. Since then the 
program has been successfully used for calculating cross sections of many 
different processes, as, e.g., all standard model (SM) processes of the form 

e + e~ -> bfj{bf 2 m, (2) 

where fx, f 2 = u e , u T , u, c and /{, f 2 = e~, r~ ,d,s j2] that are relevant 
for the associated production an decay of a top quark pair and a light Higgs 
boson at the e + e~ linear collider [3], [J]. There are 240 966 Feynman diagrams 
for the hadronic channel of ([2]) 

e + e~ —> bcsbscbb, (3) 

in the unitary gauge, assuming vanishing masses of light particles m e = m s = 
0, and neglecting the Cabibbo-Kobayashi-Maskawa (CKM) mixing. Because 
of so many Feynman diagrams, the matrix element M of process ()3]), which 
is calculated in the helicity base, is rather complicated. However, if the 
Monte Carlo (MC) summing over helicities is applied, calculating \M\ 2 is 
not a problem in practice. The main issue is to calculate the integral over 
8 x 3 — 4 = 20 dimensional phase space of process (jSJ). 



2 Phase space integration in carlomat 

The phase space integration in carlomat is performed according to the fol- 
lowing algorithm. First , final state particles {^3,^4, . . . ,p n } of process §B) 
are divided into two subsets of four momenta and q i2 each, that are de- 
fined in the relative centre of mass system (r.c.m.s.), + q^ = 0. This is 
done in a way that depends on the topology of a diagram. Then the identity 

is inserted in the standard parametrization of the Lorentz invariant phase 
space element 

n n j 3 

d^Lips = (27r)V 4 )(p 1+ p 2 -E^) II (2^- ^ 
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where n/ = n — 2 is the number of final state particles of process (CEJ). The 
insertion is repeated consecutively until parametrization (JSJ) is brought into 
the following form 

d 3n f- 4 Lips = (27r) 4 - 3 ^d/od/i...d/ n _ 4 ds 1 ds 2 ...ds n _4, (6) 

where invariants Si are given by 



(7) 



s . = l (ft'i+fta) = (^1+^2) , for 1 = 1,..., 71 -4 
\ (Pi + P2) 2 = s, for i = 

and dli are the two particle phase space elements 

2^3- w 

In Eq. (JS]), A(x, ?/, z) = x 2 + + ;z 2 — 2xy — — 2yz is the kinematical 
function and fli is the solid angle of momentum q il in the r.c.m.s. 

In carlomat v. 1.0, a different phase space parametrization ([6]) is gener- 
ated for each of N Feynman diagrams of process ([T|) 

fi(x) = d 3n f- 4 Lip Sl (x) i = l,...,N, (9) 

where x = (xi, ...,x 3nf _^j are uniformly distributed random arguments and 
the normalization condition 

1 

J dx 3n f~ 4 fi{x) = vo\{Lips) (10) 


is satisfied for each parametrization. Invariants Sj of (j7|) are randomly gener- 
ated within their physical limits which are deduced from the topology of the 
Feynman diagram. They are generated either according to the uniform distri- 
bution or, if necessary, mappings of the Breit-Wigner shape of the propaga- 
tors of unstable particles and ~ 1/s behaviour of the propagators of massless 
particles are performed. An option is included in the program that allows to 
turn on the mapping if the unstable particle decays into 2, 3, 4, ... on-shell 
particles. Different phase space parametrizations obtained in this way can 
be used for testing purposes. 
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All the parametrizations fi(x) of Eq. are then automatically combined 
into a single multichannel probability distribution 

JV 

f(x) = Y,a l f l (x), (11) 

i=l 

with non negative weights aj, i = 1, JV, satisfying the condition 

iv i 

$>i = 1 ^ / dx 3n ^ 4 /(^) = vol(Lzps). (12) 
i=i 

The actual MC integration is done with the random numbers generated ac- 
cording to probability distribution f(x). A large number of the Feynman 
diagrams, which is typical for multiparticle processes, results in the equal 
number of kinematical routines that are generated. The kinematical routines 
that contribute the most to the integral can be selected with the iterative 
approach described below. 

Integration in carlomat can be performed iteratively. First, the MC 
integral is calculated JV times with a rather small number of calls to the 
integrand, each time with a different phase space parametrization fi(x), all 
the parametrizations being assigned weight Oj = 1/JV. The result cjj obtained 
with i-th parametrization is used to calculate a new weight according to the 
following formula 

N 

a i = °i/ J2 °3 ( 13 ) 

that is the probability of choosing i-th parametrization in the first iteration. 
In this way channels with small weights a, are not chosen and will have zero 
weights in the next iteration. After the first iteration has been completed, 
the new weights for the second iteration are determined analogously and so 
on. After several iterations only the most important kinematical channels 
survive. However, the large number of kinematical channels for multiparticle 
processes in the beginning implies a very long compilation time. 

This has been improved in the current version of carlomat by introducing 
the following changes: 

• Commands for calculating phase space boundaries and boosts of the 
particle four momenta from the r.c.m.s, where they randomly gener- 
ated, to the c.m.s. are generated only once for all parametrizations 
corresponding to diagrams of the same topology. 
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• Kinematical routines corresponding to the diagrams of the same topol- 
ogy that contain the same mappings are discarded at the stage of code 
generation. 

In this way, the generated code is substantially shortened and reduction of 
a compilation time, typically by a factor 2 — 5 for multiparticle processes, is 
achieved. 



3 Extensions of SM 



Extensions of SM that have been implemented in the current version of 
carlomat include scalar electrodynamics and an anomalous Wtb coupling. 
The details on the latter are described below. 

The most general effective Lagrangian of the Wtb interaction contain- 
ing operators of dimension four and five that has been implemented in the 
program has the following form [5] 
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where mw is the mass of the W boson, Pl 



Id 



75 J 



are the left- and right-handed chirality projectors, a^ v 



;i + 7 5 ) 

I, V th 



and Pr = 

= | [f,rh Vtb is 

the element of the CKM matrix with the superscript * denoting complex 
conjugate, f^, ff 1 , ff 1 and i — 1,2, are form factors which can be com- 
plex in general. Other dimension five terms that are possible in Lagrangian 
( I14p for off shell W bosons vanish if the W's decay into massless fermions, 
which is a well satisfied approximation for fermions lighter than the 6-quark. 
The subroutines necessary for calculation of the helicity amplitudes involving 
the right- and left-handed tensor form factors of Lagrangian ( fl4l) have been 
written and thoroughly tested. 

The lowest order SM Lagrangian of the Wtb interaction is reproduced by 
setting 



ft = ft 



/* = f * = ft = f* = f 2 R = f. 







(15) 
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in ( I14p . If CP is conserved then the following relationships hold 

fT = f?, fi* = fi ^d /?' = #, M* = A R (16) 

Thus, 4 independent form factors are left in Lagrangian ( 114ft . See [0] for the 
Feynman rules resulting from (fl4l) . 

Direct Tevatron limits, obtained by investigating two form factors at a 
time and assuming the other two at their SM values, are the following [7] 

< 0.28. (17) 

The direct LHC limits are still weaker [Sj. If CP is conserved then the 
right-handed vector form factor and tensor form factors can be indirectly 
constrained from the CLEO data on b — > S7 branching fraction [9] and from 
other rare B decays [10] . However, there is still some room left within which 
the anomalous form factors, in particular the tensor ones, can be varied. 

The new version of carlomat may be useful when looking for new physics 
effects in the processes of the top quark production both at hadron-hadron 
and e + e~ collisions. In the context of the latter, a feature of the program that 
may prove itself particularly useful is that full information about helicities 
of the external particles can be easily retrieved by switching off the MC 
summing over polarization states of one or more particles. 

Let us illustrate the usefulness of the automatic approach to the imple- 
mentation of the anomalous Wtb coupling by looking closer at the top quark 
pair production in hadron-hadron collisions. Main processes of ti production 
at hadron colliders are 

qq ->■ ti, gg -> it. (18) 

The quark- ant iquark annihilation dominates at Tevatron while the gluon- 
gluon fusion dominates at LHC. Taking into account decays (t — > bW — > 
bff) results in processes with 6 particles in the final state, as e.g. 

uu —> bu^ji^ bfi'u^, (19) 
uu —> bu^fi + bdu, (20) 
uu — > bud bdu (21) 

with 559, 718 and 6134 Feynman diagrams in the lowest order, respectively 
(unitary gauge, m u = rrid = m s = m e = = 0, no CKM mixing). Ex- 
amples of the Feynman diagrams of process (|2"U|) are shown in Fig. [TJ Let 
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< 1.01, 



< 0.23, 
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Figure 1: Examples of the lowest order Feynman diagrams of process (120p . 
Black blobs indicate the Wtb coupling. 

us note that the Wtb coupling that is indicated by a black blob enters twice 
both in the ti production signal diagrams depicted in Fig. [TJi and [Tb and the 
single top production diagram of Fig. [Tb. Obviously it is not present in the 
off resonance background diagrams an example of which is shown in Fig. [TH. 
Needles to say that the by hand implementation of the anomalous Wtb cou- 
pling in the matrix element of any of processes ( !T9|) -(!20l) would have been 
rather tedious a task. With the current version of carlomat the cross sec- 
tion of these processes can be computed automatically with any choice of the 
anomalous form factors of (fill) , as it was done in [TT] , where an influence of 
the anomalous Wtb coupling on forward-backward asymmetry of top quark 
pair production at the Tevatron was investigated taking into account decays 
of the top quarks to 6 fermion final states containing one charged lepton. 

4 Summary and Outlook 

New developments in carlomat have been presented. They include improve- 
ments of the generation of phase space integration routines, implementation 
of scalar electrodynamics and the anomalous Wtb coupling including oper- 
ators of dimension up to five and many other improvements, as e.g., size 
reduction of the colour matrix, that have not been discussed in the present 
lecture. Some minor bugs in the program have been corrected too. A new 
version of the program will be released, hopefully soon. 
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